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Abstract 


The  near  field  region  of  the  turbulent  wake  behind  a  stationary 
(u)*0  RPM)  and  rotating  (u>«500  RPM)  two  dimensional  circular  cylinder  is 
documented.  U  component  mean  velocities,  turbulent  intensities  and 
velocity  autocorrelation  functions  are  measured  at  various  downstream 
locations.  From  the  autocorrelation  functions,  Taylor  time  microscales, 
longitudinal  integral  time  scales  and  energy  dissipation  rates  are  com¬ 
puted.  The  free  stream  velocity  and  Reynolds  number,  based  on  cylinder 
diameter,  are  respectively  5.9  m/sec  and  approximately  50,000.  The 
turbulent  measurements  are  made  with  a  laser  velocimeter  using  a  fre¬ 
quency  shifting  phase  modulator  and  the  photon  correlation  processing 
scheme.  Naturally  occurring  contaminant  is  augmented  with  vaporized 
kerosene  to  increase  the  magnitude  of  scattered  light.  Results  are 
presented  for  the  effects  of  varying  angular  speed,  frequency  shifting, 
particle  concentrations  and  detector  optics  pinhole  size  on  the  turbu¬ 


lent  quantities. 


I  INTRODUCTION 


Background 

Several  studies  have  been  completed  at  the  Air  Force  Institute  of 
Technology  utilizing  a  Laser  Doppler  Velocimeter  system  in  conjunction 
with  the  photon  correlating  scheme  (References  1,  2,  3  and  A).  A 
common  factor  of  all  of  these  studies  has  been  the  use  of  naturally 
occurring  flow  contaminant  to  obtain  laser  beam  light  scattering.  The 
present  study  was  conceived  as  a  further  demonstration  of  the  LDV's 
capabilities,  by  studying  the  flow  about  a  two  dimensional  circular 
cylinder  with  naturally  occurring  contaminant  augmented  with  flow 
seeding.  The  facility  chosen  for  this  investigation  already  existed 
in  the  form  of  the  Air  Force  Institute  of  Technology  two  dimensional 
smoke  tunnel.  This  facility  possesses  velocity  ranges  within  the 
capability  of  the  LDV  and  is  ideally  suited  for  marker  particle  seeding 
because  of  its  inherent  smoke  generating  system. 

The  choice  of  a  two  dimensional  circular  cylinder,  as  the  model 
for  study,  was  made  because  of  the  diversity  of  flow  regions  that 
comprise  its  flowfield  (Figures  1  and  2) .  The  turbulent  near  wake  of 
this  investigation  can  be  described  as  incompressible,  unsteady  and 
initial  condition  dependent. 

Available  theory  for  two  dimensional  wakes  is  restricted  to  the 
far  field  region,  where  self  preservation  (the  only  change  in  the  flow 
is  in  the  length  and  velocity  scales  used  to  describe  it)  and  Reynolds 
number  similarity  exist.  The  far  field  wake  does  not  begin  for  ten 
diameters  downstream  of  the  cylinder. 
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Experimental  studies  of  the  two  dimensional  wake  of  a  circular 
cylinder  are  many,  but  their  comparison  the  present  study  is  very 
restricted.  Most  work  appears  to  have  been  centered  about  very  low 
Reynolds  number  (<200)  flows.  The  wake  structure  at  these  Reynolds' 
numbers  (based  on  cylinder  diameter)  bears  little  resemblence  to  the 
present  studies  wake  structure  (Reynolds  number  based  on  cylinder 
diameter  of  50,000). 

The  far  field  wake  is  treated  to  the  practical  exclusion  of  the 
treatment  of  the  near  wake.  Measurements  extend  to  only  five  cylinder 
diameters  downstream  in  this  investigation.  The  far  wake  theoretically 
does  not  commence  until  ten  cylinder  diameters  downstream. 

In  many  cases  the  cylinder  does  not  span  the  entire  width  of  the 
test  section,  resulting  in  three  dimensional  end  effects.  In  essence 
few  empirical  results  are  available  for  exact  comparison  with  the 
measur emen  ts  taken . 

Objective 

The  main  objective  of  this  investigation  is  to  determine  the 
validity  of  using  the  AFIT  LDV  system  in  an  artifically  seeded  flow- 
field.  The  characterization  of  the  two  dimensional  circular  cylinder 
flowfleld  by  measurements  of  mean  velocity,  turbulent  intensity,  in¬ 
tegral  time  scales,  Taylor  time  microscales  and  kinetic  energy  dissipa¬ 
tion  rates,  is  secondary  in  nature,  but  intimately  tied  into  the  main 
objective.  Determination  of  the  nature  of  the  flowfield  behavior  is 
necessary  for  an  accurate  assessment  of  the  effects  of  parametric 
studies. 

Approach 

The  approach  taken  to  meet  the  main  objective  of  this  study  is  to 
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conduct  several  parametric  studies  within  the  different  flow  regions 
behind  the  two  dimensional  circular  cylinder.  The  study  is  conducted 
at  one  free  stream  velocity  (5.9  meters/second).  The  parameters  selected 
for  variation  are  marker  particle  concentration,  aperture  size,  frequency 
shift  of  the  phase  modulator  and  to  a  lesser  extent  rotation  rate  of  the 
cylinder.  The  flowfield  regions  chosen  for  probing  include  the  potential 
flow  region  above  the  cylinder,  the  boundary  layer  near  the  surface  of 
the  cylinder,  the  recirculating  region  behind  the  cylinder  and  the  near 
wake,  downstream  of  the  recirculating  region.  An  evaluation  of  the 
compatibility  of  the  LDV  in  conjunction  with  the  seeded  flow  system  will 
be  based  on  the  effect  of  parametric  variation  on  local  mean  velocity 
and  turbulent  intensity  measurements  and  how  they  agree  with  expected 
results. 

This  investigation  is  herein  documented  with  discussions  of  the 
experimental  apparatus  and  procedure  utilized,  the  results  of  the  measure¬ 
ments,  the  conclusions  drawn  from  these  measurements,  recommendations  for 
follow-on  investigations  and  appendices  of  detailed  data  reduction  methods. 
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II  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

Flow  System 

The  facility  used  for  this  investigation  is  the  Air  Force  Institute 
of  Technology  two  dimensional  smoke  tunnel  (Figure  3) •  This  facility 
employs  an  open  return  system  of  flow,  capable  of  subsonic  incompressible 
velocities  up  to  23  meters/second,  using  two  diffuser  isolated  1.5  horse¬ 
power  motor  driven  fans.  Turbulent  measurements  are  taken  at  a  nominal 
frees tream  velocity  of  5.9  m/sec  and  Reynolds  number  based  on  cylinder 
diameter  (11.75  cm)  of  approximately  50,000.  The  velocity  was  monitored 
at  first  by  a  Prandtl  type  pitot  probe  and  a  micromanometer.  These  were 
later  discarded  in  favor  of  the  laser  velocimeter,  as  a  means  of  setting 
and  checking  the  freestream  velocity. 

The  removable  front  test  section  measures  1.5  m  in  length,  1.0  m  in 
height  and  0.07  m  in  width.  The  back  wall  of  the  test  section  is  of 
laminated  plate  glass,  whereas  the  front  wall  is  of  0.0097  m  thick  plexi¬ 
glass.  This  window  arrangement  is  acceptable  in  light  of  the  fact  that 
the  laser  velocimeter  is  operated  in  the  backscatter  mode.  The  test 
section  is  noted  as  having  a  downhill  gradient  of  0.05  m  in  1.5  ra. 

The  flow  marker  particles  are  introduced  into  the  flow  system  in 
thin  streamtubes  by  a  stack  and  injector  apparatus,  positioned  in  the 
tunnel  contraction  region.  The  stack  is  of  airfoil  shape  with  sixty 
five  0.60  cm  inside  diameter  injector  tubes  issuing  from  it.  The  marker 
particles  are  generated  in  a  process  where,  two  900  watt  Inconel  heaters 
boll  kerosene  fuel  at  605  degrees  Kelvin,  creating  vapor  particles 
that  are  then  mixed  under  turbulent  conditions  with  cool  air  to 
produce  a  dense  white  nontoxic  and  noncorrosive  smoke.  Water 
vapor  is  condensed  out  of  the  smoke  in  a  condensing  chamber  prior 
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to  entrance  into  the  stack/injector  apparatus,  thereby  eliminating  water 
condensation  in  the  injector  tubes  and  the  test  section. 

A  great  deal  of  prior  effort  was  spent  in  reducing  turbulence 
levels  in  this  flow  system  (Reference  5).  These  efforts  included  in¬ 
stallation  of  a  0.2  m  radius  bell  mouth  to  the  tunnel  inlet,  in  an 
attempt  to  correct  for  a  low  contraction  ratio  of  11.5  to  1.0.  Immedi¬ 
ately  downstream  of  this  location  a  0.076  m  thick  section  of  honeycomb, 
with  a  cell  length  to  diameter  ratio  of  8.0,  was  installed  to  reduce 
large  scale  turbulent  structures.  In  addition  a  series  of  screens 
(progressively  finer  downstream)  were  mounted  prior  to  the  stack/ 
injector  location.  From  measurements  taken  in  the  freestream,  it  would 
appear  that  turbulent  intensities  approach  a  value  of  0.10  in  the  smoke 
stream tubes.  Such  high  values  must  be  attributed  to  the  shape  of  the 
stack/ injector  system,  its  location  and  the  process  of  issuing  a 
secondary  flow  into  the  mainstream  through  the  injectors  (Reference  5), 
Circular  Cylinder 

Two  circular  cylinders,  one  capable  of  rotation  and  one  not,  were 
constructed  that  mounted  to  the  test  section  back  wall  in  the  approxi¬ 
mate  center  (Figure  3).  Both  cylinders  are  of  0.059  m  radius  and  0.07 
m  thickness.  This  represents  a  test  section  area  blockage  ratio  of 
12%.  One  cylinder  is  constructed  of  polished  wood  with  a  static  pres¬ 
sure  orifice  located  at  midspan.  A  positive  seal  with  the  test  section 
sidewalls  is  achieved  by  capping  each  end  of  the  cylinder  with  felt, 
thereby  assuring  two  dimensional  flow  about  the  cylinder.  A  cylinder 
surface  pressure  survey  was  conducted  (Figure  4)  at  two  values  of  Upg 
and  is  compared  with  invlscld  theory  and  experimental  measurements  for 
8ubcritlcal  flows  (Reference  6). 
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The  cylinder  capable  of  rotation  (Figure  5)  is  constructed  of 
aluminum,  with  a  smooth  surface  painted  flat  black  to  minimize  scatter¬ 
ing  of  light  from  its  surface.  Two  endplates  of  plexiglass,  one  fixed 
and  one  free  floating  are  mounted  to  either  end.  With  the  cylinder  in¬ 
stalled  in  the  test  section  a  gap  of  approximately  0.051  mm  exists  be¬ 
tween  the  plexiglass  endplates  and  aluminum  main  body.  The  possibility 
of  three  dimensional  flow  existing  at  the  gap  in  the  cylinder  ends  is 
considered  negligible.  Rotation  of  the  cylinder  is  accomplished  by 
mounting  a  small  variable  speed  motor  to  a  shaft  extending  from  the 
cylinder  main  body  through  the  back  wall.  Rotations  of  up  to  u>=800  RPM 
are  possible  without  any  significant  vibrations  being  introduced  into 
the  test  system.  Rotation  rates  are  constant  to  within  +  5  RPM. 

Laser  Velocimeter 

The  laser  used  for  all  measurements  is  a  Helium-Neon  Laser  of  15 

milliwatts  intensity  at  6328.0  x  10  ^  m,  plus  associated  power  unit. 

2 

The  laser  beam  diameter  is  1.1  mm  at  the  1/e  points.  The  transmitting 
optics  consisted  of  a  transmitter  beamsplitter  and  polarization  unit 
mounted  to  the  laser  head,  a  frequency  shifting  electro-optic  phase 
modulating  crystal,  two  front  surface  silvered  plane  mirrors  mounted  at 
45  degrees  to  the  horizontal  and  a  convex  focussing  lens  of  100  cm  focal 
length.  The  beamsplitter  takes  the  incident  beam  from  the  laser  and 
divides  it  into  two  equal  intensity,  1.1  mm  diameter  beams.  The  prisms 
used  in  this  process  are  anti-reflection  coated  at  6328.0  x  10  ^  m. 

Two  adjustments  to  the  laser  beams  are  available.  First,  there  is  an 
adjustment  for  a  (the  half  angle  between  the  laser  beams  at  the  control 
volume).  For  this  investigation  the  value  of  a  is  equal  to  a  constant 
0.6  degrees.  Second,  there  is  an  adjustment  for  beam  separation.  A 
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constant  value  of  £*20  mm  is  maintained. 


The  frequency  shifting  phase  modulator  is  required  in  the  optical 
train  to  eliminate  any  flow  direction  sense  ambiguities  and  provide 
measurements  in  regions  of  high  turbulent  intensity.  In  principle  a 
uniform  shifting  of  the  fringes  in  the  control  volume  (point  of  laser 
beam  intersection)  is  possible  with  an  application  of  a  sawtooth  volt¬ 
age  to  the  phase  modulator’s  two  electro-optical  crystals.  A  resulting 
increase  or  decrease  in  the  doppler  frequency  of  the  flow  enables  the 
flow  direction  to  be  determined.  A  drive  unit  is  required  for  the 
phase  modulator  unit  as  well  as  a  universal  counter  for  accurate  de¬ 
termination  of  the  doppler  frequency  shift  applied. 

Two  front  surface  silvered  mirrors  are  required  for  transmission 
of  the  laser  beams  along  the  optical  path.  Both  mirror  surfaces  are 
masked  to  eliminate  stray  light  from  being  transmitted  through  the 
system.  The  convex  focussing  lens  is  made  of  high  quality  optical 
glass. 

The  receiving  optics  consists  of  a  collecting  lens/aperture  and 
photo  multiplier  tube.  The  receiver  photon  detection  unit  employs  a 
200  mm  telephoto  lens  as  a  spatial  filter  by  being  focused,  with  an 
external  reflex  viewer,  onto  the  fringe  control  volume.  The  scattered 
light,  from  particles  passing  the  alternate  constructive  and  destructive 
interference  fringes  of  the  control  volume,  are  foe  ised  onto  a  400  pm 
diameter  pinhole  prior  to  detection  by  the  photo  multiplier  tube.  The 
400  ym  diameter  pinhole  provides  the  best  signal  to  noise  ratio  for  the 
current  experimental  setup.  The  associated  photo  multiplier  tube  power 
unit  supplies  a  constant  1850  volts  to  the  cathode  of  the  photomultiplier 
tube. 
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The  laser,  transmitting  optics  and  receiving  optics  are  all  mounted 
to  a  wood  and  metal  optical  bench  (Figures  6  and  7).  Flowfield  longi¬ 
tudinal  traverses  are  accomplished  by  translation  of  the  entire  bench. 
Lateral  traverses  are  conducted  with  an  aluminum  secondary  table  capable 
of  uniform  translation  through  a  chain  and  sprocket  mechanism.  The 
design  of  the  secondary  table  limited  lateral  traverses  in  the  -Y 
direction  to  Y/2Ro®-1.0.  A  traverse  across  the  test  section  in  the 
Z  direction  was  conducted  early  in  the  investigation  to  determine  the 
extent  of  twp  dimensional  flow  across  the  width.  The  resultant  velocity 
profile  indicated  a  uniform  velocity  2.0  cm  on  either  side  of  the  test 
section  centerline.  All  measurements  were  conducted  at  the  centerplane 
of  the  test  section. 

The  electronic  processing  scheme  is  composed  of  a  digital  photon 
correlator  and  data  storage  unit,  associated  power  supply  and  oscilli- 
scope  for  visual  observation  of  the  autocorrelation  function  growth 
with  time.  The  correlator  possesses  a  resolution  time  of  from  50  nano¬ 
seconds  to  1  second.  Measurements  were  taken  in  the  single  clipped 
autocorrelation  mode  and  at  ar  infinite  sample  rate.  A  schematic  of 
the  electronic  signal  processing  equipment  is  shown  in  Figure  8* 


13 


VELOCIMETER  SETUP 


Ill  RESULTS  AND  DISCUSSION 


Measurements  are  made  in  the  near  wake  of  a  two  dimensional  cir¬ 
cular  cylinder  of  one  point  statistical  properties.  These  properties 
consist  of  mean  velocities,  turbulent  intensities  and  autocorrelation 
functions  in  the  X  direction.  A  sample  of  typical  data  is  shown  in 
Table  I. 


R( VOLTS,  ARBITRARY) 

T(ySECONDS) 

2*08  lst  RMAXIMA 

0.0 

2.06 

.0031 

2.02 

.0063 

1.96 

.0094 

1.90 

.0125 

1.84 

.0156 

1.80 

.0188 

1.77 

.0219 

1<76  lst  Minima 

.0250 

1-98  2nd  rmaxima 

.0531 

1.72  2nd 

.0781 

i.88  3rd  I*MAXIMA 

.1094 

1,64  3rd  ^NIMA 

.1375 

1.77  4th  Rmaxima 

.1656 

1.56  4  th 

.1969 

TABLE  I  Digital  Autocorrelation  Function  Data, 
X/2Ro— 3,13,  Y/2RO-0.09 


This  data  is  a  digital  representation  of  a  typical  autocorrelation  func¬ 
tion  (see  Figure  37).  It  consists  of  hand  recorded  values  of  R(auto- 
correlation  function  ordinate  in  arbitrary  volts)  versus  ^autocorrela¬ 
tion  function  abscissa  in  ^seconds,  i.e.  delay  time).  From  this  data 
longitudinal  mean  velocities  and  turbulent  intensities  are  computed  as 


? 

detailed  in  Appendix  A,  If  the  autocorrelation  function  is  skewed  (as 
it  was  for  this  entire  investigation)  it  must  be  corrected  for  skewness 
as  shown  in  Appendix  B.  In  addition  to  the  skewness  correction  a 
correction  for  ambiguity  noise  must  be  applied  to  the  autocorrelation 
function.  This  procedure  is  outlined  in  Appendix  C.  Integral  time 
scales  and  Taylor  time  microscales  are  also  computed  from  the  skewness 
corrected  autocorrelation  function.  These  procedures  are  treated  in 
Appendices  D  and  E  respectively.  Finally,  kinetic  energy  dissipation 
rates  are  estimated  from  the  integral  time  scales  and  Taylor  time 
microscales . 

The  effects  of  cylinder  rotation,  marker  particle  concentration, 
detector  pinhole  diameter  and  frequency  shifting  on  the  mean  velocities 
and  turbulent  intensities  are  noted. 

A  laser  doppler  velociroeter  is  used  to  obtain  all  measurements. 

Mean  Velocities 

Lateral  mean  velocity  profiles  are  presented  in  Figures  9  to  16 
for  downstream  locations  X/2Ro«0.0,  0.25,  0.50,  0.75,  1.00,  2.00,  3.00 
and  5.00,  with  the  cylinder  both  rotating  and  nonrotating.  The  profiles 
are  assumed  to  be  symmetrical  about  the  wake  centerline  for  the  station¬ 
ary  case. 

The  local  mean  velocity  is  nondimens ionali zed  by  the  free  stream 
velocity  and  the  lateral  displacement  is  nondimensionalized  by  twice 
the  cylinder  radius. 

Note  the  existence  of  a  region  of  reverse  flow  behind  the  circular 
cylinder  in  Figures  9,  10  and  11.  This  recirculating  region  begins  at 
X/2Ro«0.0  and  decreases  in  width  with  increasing  downstream  location, 
until  it  terminates  somewhere  between  X/2Ro«0.75  and  X/2Ro*1.00 
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LATERAL  MEAN  VELOCITY 


(Figure  13),  Observe  that  large  velocity  gradients  (3U/3Y)  exist  at 
the  boundaries  of  the  recirculating  region.  The  local  mean  velocities 
just  outside  of  the  recirculating  region  exhibit  values  in  excess  of 
freestream  (i.e.  U/UFS>1.0).  This  is  attributed  to  the  fact  that  the 
flow  has  experienced  an  acceleration  and  then  deceleration  in  traver¬ 
sing  around  the  cylinder.  The  deceleration  to  Up^  with  downstream 
location  (U/UF$*1.4  at  X/2Ro=0.0,  whereas  U/UFS*1.2  at  X/2Ro=0.50)  in 
the  area  of  the  recirculating  region  is  not  complete.  Effects  of 
cylinder  rotation  at  u>=500  RPM  on  mean  velocities  are  noted  as  being 
negligible  in  the  recirculating  region. 

At  further  downstream  locations  (X/2Ro=0.75  to  X/2Ro=5.00  or 

Figures  12  to  16)  note  that  the  profile  velocity  deficit  (UWA  -UT -^AT ) / 

MAX  LOCAL 

U  decreases  in  the  downstream  direction  and  the  wake  width  appears  to 
Jr  b 

be  growing.  Figure  17  addresses  the  apparent  growth  of  the  wake  half 
width.  Wake  half  width  is  defined  as: 

(U/U  )  -  (U/U  ) 

Ym=  MAX _  CENTERLINE 

2 

It  is  nondimensionalized  by  the  cylinder  radius  and  plotted  versus  the 
downstream  location  nondimensionalized  by  twice  the  cylinder  radius. 
Neglecting  the  data  between  X/2Ro=0.0  and  X/2Ro=2.00  (dashed  line)  it 
appears  that  the  wake  half  width  is  growing  at  a  rate  of  approximately 
5  degrees. 

At  the  downstream  locations  X/2Ro=2.0  and  X/2Ro=3.0  note  the 
existence  of  a  region  of  higher  velocity  flow  at  the  centerline  of  the 
wake  and  of  lower  velocity  flow  on  either  side  of  the  centerline.  This 
feature  is  more  distinct  at  X/2Ro*2.0  than  X/2Ro«3.0.  Therefore*  it 
appears  to  be  decaying  with  downstream  location.  Consideration  of 
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location  X/2Ro*5.0  (Figure  16)  reveals  that  this  is  the  case.  This 
behavior  will  be  discussed  in  more  detail  when  variation  of  centerline 
velocity  is  treated. 

Returning  again  to  the  effects  of  cylinder  rotation  on  the  mean 
velocities,  a  scan  of  the  mean  velocity  profiles  downstream  of  the 
recirculating  region  (Figures  12  to  16)  reveal  that  rotation  at  u)=500 
RPM  has  no  effect  here,  as  was  also  the  case  in  the  recirculating 
region. 

The  development  of  centerline  mean  velocities  downstream  is  shown 
in  Figure  18.  Note  the  recirculating  region  with  its  reverse  flow  and 
the  immediately  following  steep  velocity  gradient.  The  return  to  free- 
stream  velocities  (dashed  line)  at  centerline  is  rapid  (within 
X/2Ro=4.00)  as  expected  in  flow  with  a  low  freestream  velocity  of 
5.9  m/sec.  A  slight  overshoot  in  centerline  velocity  occurs  between 
X/2Ro=1.00  and  X/2Ro~3.00,  followed  by  an  asymptotic  decay  to  free¬ 
stream  values. 

One  of  the  parameters  of  interest  in  this  investigation  is  the 
effect  of  cylinder  rotation  rates  on  the  mean  velocities  in  the  cylinder 
flowfield.  The  preceeding  discussion  revealed  the  lack  of  effect  of 
rotation  on  the  cylinder  near  wake.  Further  evidence  of  this  is  pro¬ 
vided  in  Figure  19.  The  location  X/2Ro*0.0,  Y/2Ro=0.34  is  in  the 
cylinder  recirculating  region.  For  values  of  u)  from  0  to  800  RPM  no 
significant  change  in  U/U  c  is  noted.  Also  in  Figure  19  the  effect  of 

r  b 

rotation  for  location  X/2Ro*=0.75,  Y/2Ro*0.6  is  presented.  This  location 
is  considered  to  be  significant,  because  it  is  located  outside  of  the 
recirculating  region  and  in  the  region  of  vortex  formation.  The  effect 
of  rotation  is  again  insignificant. 
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FI6URE  19.  EFFECT  OF  CYLINDER  ROTATION  ON  MEAN  VELOCITIES 


31 


3 


The  effect  of  rotation  on  the  mean  velocities  outside  of  the  near 
wake  and  in  a  region  of  potential  flow  is  shown  in  Figure  20.  A  lateral 
traverse  in  the  positive  Y/2Ro  direction  is  shown  for  X/2Ro*-0.50. 
According  to  potential  flow  theory  at  9*90  degrees,  the  local  velocity 
should  increase  as  the  cylinder  surface  is  approached  radially  according 
to 


until  a  value  of  2  U  is  reached  at  A*Ro,  9*90  degrees  (the  cylinder 

F  b 

surface).  Agreement  between  the  potential  flow  theory  and  the  measured 
mean  velocities  at  a>=  0  RPM  is  reasonable.  With  u>*500  and  rotating  clock¬ 
wise  an  increase  in  local  mean  velocities  is  noted  from  the  surface 
(Y/2Ro*0.5)  out  to  Y/2Ro*1.36.  Reversing  the  direction  of  rotation  (from 
clockwise  to  counterclockwise)  produces  lower  local  mean  velocities  out 
to  Y/2Ro=1.36.  It  would  appear  then,  that  the  effects  of  rotation  on 
mean  velocities  are  most  significant  in  potential  flow  regions  close  to 
the  cylinder  surface.  Flow  acceleration  around  the  cylinder  is  again  in 
evidence  (at  X/2Ro*-0.5)  by  the  return  of  U/Upg  to  values  above  1.0  for 
lateral  positions  increasing  outward  from  the  cylinder.  The  ability  of 
the  laser  doppler  velocimeter  to  define  the  cylinder  boundary  layer  is 

shown  by  the  rapid  decrease  in  mean  velocity  from  U/IL*1.8  to  U/U-*0.0 

r  b  r  o 

in  a  AY/2Ro  of  0.045.  This  corresponds  to  a  boundary  layer  thickness  of 
slightly  over  one  half  of  a  centimeter. 

A  second  parameter  of  interest  was  the  effect  of  marker  particle 
concentrations  on  mean  velocity  measurements*  High  levels  of  marker 
particle  concentration  can  result  in  two  phase  flows,  creating  a  flow- 
field  different  from  the  desired  one.  Marker  particle  concentrations 
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X/2Ro=  -0.5 


-Aft* 


U/UFS 

*w/u 

0.99 

0.15 

1.05 

— 

0.94 

0.14 

0.94 

0.15 

were  not  qualitatively  determined  in  this  investigation.  It  was  deemed 
sufficient  to  visually  observe  the  intensity  of  scattered  laser  beam 
light  coming  from  the  control  volume  and  assume  that  higher  intensities 
indicate  higher  concentrations.  Actually,  the  relation  between  visual 
observation  and  intensity  is  a  logarithmic  one.  The  calculated  mean 
velocities  and  turbulent  intensities  (Appendix  A)  based  upon  this  con¬ 
centration  study  are  presented  in  Table  II. 


MARKER  PARTICLE 
CONCENTRATION  LEVEL 


1  (low) 

2 

3 

4  (high) 


TABLE  II  Effect  of  Marker  Particle  Concentration  on 
Mean  Velocities  and  Turbulent  Intensities 
at  X/2Ro=5.00,  Y/2Ro=0.0. 

Any  deviations  are  within  the  experimental  error  (approximately  5 
percent)  and  the  concentration  levels  in  use  are  assumed  to  produce  no 
two  phase  flow  effects. 

A  factor  that  greatly  affects  the  reliability  of  laser  doppler 
velocimeter  measurements  is  the  size  of  the  pinhole  diameter  in  the 
detector  system.  Three  pinhole  diameters  are  available  for  study  (100 
pm,  200  pm  and  400  pm).  The  effect,  on  mean  velocities  in  two  differing 
flowfield  regions,  of  pinhole  diameter  was  studied.  Table  III  a  pre¬ 
sents  results  taken  in  the  cylinder  recirculating  region  at  X/2Ro*0.50, 
Y/2Ro*0.0.  U/Upg  values  are  unchanged  with  varying  pinhole  diameter. 
Tables  III  b  and  III  c  present  results  in  the  flowfield  region  upstream 
of  the  cylinder  (free  stream).  Any  variations  detected  in  U/U„c 
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approach  the  value  of  the  experimental  error  (approximately  5%), 

Finally,  little  information  exists  on  the  effect  of  frequency 
shifting  on  turbulence  measurements  (Reference  7).  In  this  investiga¬ 
tion  regions  of  very  high  turbulence  (greater  than  30%)  were  encounter¬ 
ed,  which  necessitated  the  use  of  the  frequency  shifting  crystal  os¬ 
cillator.  A  study  was  made  to  determine  the  effect  on  mean  velocities 
of  this  shifting  and  the  results  are  presented  in  Figure  21.  The 
cylinder  recirculating  region  and  several  locations  upstream  of  the 
cylinder  were  probed.  Mean  velocities  are  plotted  versus  the  amount  of 
frequency  shift  in  KHz.  Both  regions  exhibit  a  nonlinear  trend  with 
increasing  shift.  Frequency  shifts  employed  in  the  freestream  region 
ahead  of  the  cylinder  of  up  to  20  KHz  appear  to  be  correct,  since  the 
values  of  U/UT7C  nonshifted  and  1I/U-,-  shifted  are  both  close  to  1.0.  In 

r  b  ro 

the  cylinder  recirculation  region  shifts  of  from  200  KHz  to  800  KHz 
appear  to  produce  results  that  are  more  in  agreement  with  the  low  nega¬ 
tive  velocities  expected.  The  vast  majority  of  measurements  taken  in 
this  investigation,  when  the  frequency  shifting  was  required,  were  at 
200  KHz. 

Turbulent  Intensities 

Turbulent  intensity  is  defined  as  the  ratio  of  the  root  mean  square 
of  the  turbulent  velocity  fluctuations  to  a  reference  mean  velocity. 

In  this  investigation  the  reference  mean  velocity  is  the  free  stream 
velocity .  Therefore , 

Turbulent  intensity  *  urms^FS 

where:  urms  “  /  U'2  dt^/2 

o 
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All  turbulent  intensity  measurements  have  been  corrected  for 
ambiguity  noise.  Ambiguity  noise  is  defined  as  an  apparent  turbulence 
resulting  from  the  finite  transit  times  of  marking  particles  as  they 
pass  through  the  control  volume  fringes  in  a  random  manner.  It  is  a 
direct  result  of  the  finite  size  of  the  laser  measuring  volume 
(Appendix  C) .  This  correction,  which  is  a  constant  for  the  investiga¬ 
tion  optical  setup,  is  subtracted  from  the  computed  turbulent  intensi¬ 
ties.  Its  value  was  determined  to  be  0.039  or  3.9%.  Accuracy  of  turbu¬ 
lent  intensity  measurements  is  approximately  5%,  since  it  is  computed 

38  URMS/UL0CAL  and  UL0CAL  iS  accurate  t0  within  5%* 

Lateral  profiles  of  turbulent  intensity  are  presented  in  Figures 

22  to  29  for  downstream  locations  X/2Ro=0.0,  0.25,  0.50,  0.75,  1.00, 

2.00,  3.00  and  5.00,  with  the  cylinder  both  rotating  and  nonrotating. 

The  lateral  displacement  from  the  flow  centerline  is  again  nondimens ion- 

alized  by  twice  the  cylinder  radius.  At  X/2Ro=0.0  (Figure  22)  note  the 

higher  turbulent  intensity  at  the  centerline  and  the  inclusion  of  the 

majority  of  data  between  turbulent  intensity  values  of  0.10  and  0.30 

(10%  and  30%  respectively).  Values  of  turbulent  intensity  at  locations 

upstream  of  the  cylinder  were  typically  on  the  order  of  0.10  (based  on 

an  average  of  several  measurements).  Progressing  downstream  in  the 

cylinder  recirculation  region  to  X/2Ro*0.25,  0.50  and  0.75  extremely 

high  turbulent  intensities  are  encountered.  Turbulent  intensity  values 

computed  to  be  in  excess  of  100%  are  not  plotted.  Instead  a  perceived 

trend  (for  the  u>*0  RPM  case)  to  higher  values  is  indicated  by  dashed 

lines.  The  trends  are  symmetrical  about  the  centerline.  Upon  comparison 

with  mean  velocity  profiles  at  the  same  X/2Ro  locations  (Figures  9  to 

16),  it  beccmes  apparent  that,  the  lateral  positions  (Y/2Ro)  of  the  high 


G  6J=0  RPM 
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mean  velocity  gradients  match  the  lateral  positions  of  the  high  turbu¬ 
lent  intensities.  Consider  that  the  transfer  term  in  the  equation  for 
the  kinetic  energy  of  mean  motion  (Reference  8) 

p  U"  V"  3U 
3Y 

represents  the  exchange  of  energy  from  the  mean  flowfield  to  the  turbu¬ 
lent  flowfield  by  the  interaction  of  Reynolds1  stresses  (p  U'  V")  with 
the  mean  velocity  gradient  (3U/3Y).  In  this  investigation  the  Reynolds' 
stresses  are  unknown,  but  the  mean  velocity  gradient  is  known.  The  fact 
that  the  mean  velocity  gradient  is  large  indicates  that  a  substantial 
transfer  of  energy  is  taking  place  from  the  mean  flowfield  to  the  turbu¬ 
lent  flowfield.  Therefore,  large  turbulent  intensities  are  to  be  ex¬ 
pected  in  regions  of  high  mean  velocity  gradients. 

The  scatter  in  the  turbulent  intensity  data  in  the  recirculating 
region  for  a  cylinder  rotation  rate  of  u>=500  RPM  is  difficult  to  define 
due  to  the  large  gradients.  In  regions  of  lower  gradients  the  scatter 
appears  to  be  on  the  order  of  10%.  The  only  conclusion  that  can  be 
reasonably  drawn,  as  regards  the  effects  of  rotation,  is  that  the  ex¬ 
cessive  turbulent  intensities  (>100%)  observed  for  the  u>*0  RPM  case  are 
not  observed  in  the  u>=500  RPM  case. 

Moving  downstream  from  X/2Ro*1.00  to  X/2Ro*5.00  (Figures  26  to  29), 
it  is  still  the  case  that  regions  of  highest  turbulent  intensity  corres¬ 
pond  to  regions  of  highest  mean  velocity  gradients.  The  peak  magnitudes 
observed  in  the  recirculating  region  are  not  repeated  here,  but  are 
still  significantly  above  free  stream  values  (U_¥J_/Ut«0. 10) .  By  the 
time  location  X/2Ro*5.00  is  reached  the  turbulent  intensities  have 
returned  to  near  free  stream  values.  A  check  of  the  mean  velocity 
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gradients  at  X2Ro=5.00  (Figure  29)  reveals  a  substantial  reduction  in 
their  magnitude  signalling  a  reduction  in  peak  turbulent  intensities 
as  well. 

The  effects  of  cylinder  rotation  downstream  of  the  cylinder  re¬ 
circulating  region  appear  to  be  confined  to  a  reduction  in  peak  magni¬ 
tudes  as  was  also  the  case  in  the  recirculating  region. 

The  development  of  centerline  turbulent  intensities  downstream  is 
shown  in  Figure  30.  An  additional  effect  of  cylinder  rotation  now  comes 
to  light.  In  the  recirculation  region  and  up  to  X/2Ro^2.00  a  shift  in 
the  location  of  peak  turbulent  intensity  occurs  from  X/2Ro=0.50  to 
X/2Ro=0.75.  In  addition  there  is  a  broadening  of  the  curve  between 
X/2Ro=0.00  and  X/2Ro=2.00  indicating  a  growth  of  the  region  of  higher 
turbulent  intensities  downstream.  A  rapid  return  toward  free  stream 
values  is  noted  beginning  at  X/2Ro=2.00.  The  effects  of  rotation  from 
X/2Ro=2.00  on  downstream  are  negligible. 

As  stated  previously  one  parameter  of  interest  was  the  effect  of 
varying  cylinder  rotation  rates  (u>)  on  turbulent  measurements.  The 
effect  of  rotation  on  turbulent  intensities  at  two  differing  locations 
is  shown  in  Figure  31.  At  X/2Ro*0.0,  Y/2Ro-0.34  (the  recirculating 
region)  an  increase  in  clockwise  U)  produces  a  decrease  in  turbulent 
Intensities  up  to  u>*800  RPM.  This  agrees  with  the  early  observations 
of  Prandtl  (Reference  9) . 

This  result  corresponds  well  with  that  observed  in  Figure  22  for 
X/2RO-0.0  and  u>«500  RPM.  At  X/2Ro«0.75,  Y/2Ro«0.6  (region  of  shed 
vortex  formation)  the  effect  of  increasing  u)  is  again  one  of  initial 
decrease  in  turbulent  intensity  accompanied  by  an  increase  toward 
U)«0  RPM  values  of  turbulent  intensity  at  U)*700  RPM.  This  effect  has 
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been  observed  in  other  work  with  rotating  shapes  (Reference  10)*  It  is 
then  apparent  that,  in  the  near  wake  region  of  the  circular  cylinder, 
rotation  has  insignificant  effects  upon  mean  velocities,  but  has  pro¬ 
nounced  effect  upon  turbulent  intensities.  Returning  to  the  potential 
flow  region  above  the  cylinder  (9*90  degrees),  the  lateral  variation  of 
turbulent  intensity  is  presented  in  Figure  32.  Effects  of  rotation  lie 
within  experimental  error  and  are  therefore  ignored.  Rotational  effects 
outside  of  the  boundary  layer  are  insignificant.  Note  that  free  stream 
values  of  turbulent  intensity  are  indicated  by  the  dashed  line. 

The  effect  of  marker  particle  concentrations  on  turbulent  intensi¬ 
ties  at  X/2Ro=*5.00,  Y/2Ro=0.0  are  presented  in  Table  II.  Turbulent  in¬ 
tensity  values  at  concentration  level  2  are  undefined.  As  with  the  mean 
velocities,  no  effects  of  concentration  levels  are  noted  on  turbulent 
intensities. 

Pinhole  diameter  sizes  have  a  pronounced  effect  on  turbulent  in¬ 
tensity  measurements  as  shown  in  Table  III.  In  the  recirculating  region 
(X/2Ro*0.50,  Y/2Ro*0.0)  the  smaller  pinhole  diameters  (100  ym  and  200 
ym)  give  values  of  turbulent  intensity  in  excess  of  100%,  whereas,  in 
the  freestream  locations  X/2Ro=-2.65,  Y/2Ro=0.13  and  X/2Ro=-2.13, 
Y/2Ro*-0.09  they  give  smaller  turbulent  intensities  than  the  400  ym 
diameter  pinhole.  With  smaller  pinhole  diameters,  fewer  marker  parti¬ 
cles  are  sensed  as  they  traverse  the  control  volume  fringes.  This 
results  in  a  low  signal  to  noise  ratio  entering  the  photon  correlator. 

To  the  correlator  the  signal  appears  as  white  noise.  An  autocorrelation 
of  white  noise  appears  as  a  spike  at  T*0.0  (zero  time  of  averaging)  re¬ 
sulting  in  large  turbulent  intensity  values.  In  the  recirculation  re¬ 
gion,  where  fewer  particles  would  be  expected  to  traverse  the  fringe 
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pattern,  high  turbulent  intensities  would  be  computed.  This  phenomena 
does  not  explain  the  decrease  in  turbulent  intensity  with  decrease  in 
pinhole  diameter  observed  in  the  frees tr earn.  Further  study  of  pinhole 
diameter  effects  in  the  frees tream  is  required  to  more  accurately  define 
the  causes  of  this  anomaly.  Other  measurements  of  turbulent  intensities 
in  the  freestream,  with  the  400  \rm  diameter  pinhole  resulted  in  values 
of  of  0.13,  0.11,  0.03  and  -0.02.  The  results  appear  to  be 

dependent  on  the  positioning  of  the  control  volume  in  the  streamtubes 
issuing  from  the  stack/injector  system.  An  average  value  of 
of  0.10  was  assumed  to  be  representative  of  the  freestream. 

Finally,  frequency  shift  effects  on  turbulent  intensities  are 
shown  in  Figure  33.  For  the  freestream  measurements,  which  normally  did 
not  require  frequency  shifting,  any  applied  shift  of  over  20  KHz  gives 
results  inconsistent  with  the  efforts  to  reduce  the  tunnel  turbulence 
levels  (Reference  5).  In  the  recirculation  region,  shifts  below  200  KHz 
produced  autocorrelation  functions  whose  maxima  and  minima  points  were 
undetectable.  Shifts  above  200  KHz  produce  varying  results.  Since 
200  KHz  appeared  to  produce  realistic  mean  velocity  measurements,  it  is 
assumed  that  turbulent  intensity  measurements  taken  at  a  shift  of  200 
KHz  are  also  acceptable.  Further  experimentation  and  analysis  appears 
warranted  in  order  to  better  define  the  effects  of  frequency  shifting  on 
turbulent  measurements. 

Integral  Time  Scales 

Longitudinal  integral  time  scales  are  defined  as  follows: 

T 

C 

TJS-  J  IS  U"(T)U'(T')dT" 
o 
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r  IS 

TIS=  J  R(T')dT" 


where:  R(T^)  =  autocorrelation  function 

T^,  =  time  to  first  zero  crossing  of  autocorrelation 

IS  function 

The  integral  in  this  definition  represents  the  area  under  the  auto¬ 
correlation  function  to  its  first  zero  crossing.  This  standard  method 
of  determining  integral  scales  is  departed  from  in  this  investigation. 
Details  of  the  actual  method  employed  are  given  in  Appendix  D. 

Time  scales  are  presented  instead  of  length  scales  because  Taylor* s 


hypothesis 


u  ix  (*>  -  h 


relating  length  scales  to  time  scales  cannot  be  realistically  utilized 
in  this  investigation.  The  local  mean  velocity  (U)  is  not  constant  and 
the  statistical  property  Q  is  not  stationary  for  statistical  equilibrium 
to  exist. 

Integral  scales  represent  the  largest  turbulent  scale  in  the  flow- 
field.  The  development  of  integral  time  scales  with  downstream  location, 
for  various  lateral  positions  (Y/2Ro-0.0,  1.02,  2.04  and  -0.94)  are  pre¬ 
sented  in  Figure  34  and  a  complete  compilation  is  provided  in  Table  IV. 
The  method  of  determining  the  time  scales  (Appendix  D)  is  an  approximate 
technique,  subject  to  large  errors,  as  seen  by  the  scatter  in  Figure  34. 
Ignoring  points  of  large  deviation,  it  is  possible  to  construct  a  linear 
curve  fit  through  all  available  data  points.  The  resultant  curve  in¬ 
dicates  a  slow  rate  of  increase  in  integral  time  scales  with  downstream 
location.  Further  observations  concerning  the  data  are  not  warranted 


*****  .  ^ 
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Y/2Ro=0,  <*)  =0  RPM 
Y/2Ro=0,  W=500  RPM 
Y/2Ro=l,  U)  =0  RPM 


because  of  the  approximate  nature  of  the  integral  time  scale  method  of 
determination . 

Taylor  Time  Microscales 

Taylor  microscales  are  defined  as  the  ratio  of  the  dissipation  of 
energy  to  the  total  amount  of  energy  in  the  flow.  They  lie  between  the 
largest  scales  (integral  scales)  and  the  smallest  scales  (Kolmogorov 
scales).  Their  experimental  determination  is  outlined  in  Appendix  E. 

The  centerline  variation  of  Taylor  time  microscales  with  downstream 
location  is  addressed  in  Figure  35.  A  complete  tabulation  of  variations 
with  downstream  location,  for  several  lateral  positions  is  given  in 
Table  V.  Returning  to  Figure  35,  if  large  deviation  points  are  ignored, 
a  curve  can  be  fitted  to  the  data  that  indicates  a  slight  increase  in 
Taylor  time  microscales  downstream.  The  growth  rate  of  this  curve  is 
approximately  equal  to  that  of  the  integral  time  scale.  Cylinder  rota¬ 
tion  has  no  appreciable  effect  on  this  curve. 

Kinetic  Energy  Dissipation  Rates 

The  kinetic  energy  dissipation  rate  is  proportional  to  the  Taylor 
velocity  microscale  cubed,  over  the  integral  length  scale  (Reference  11): 


where : 


TIS  xu 


J 


In  order  to  calculate  the  kinetic  energy  dissipation  rate,  it  is 
necessary  to  assume  that  the  Taylor  hypothesis  holds.  With  this  assump¬ 
tion  made  and  recalling  the  approximate  nature  of  the  integral  time  scale 
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FIGURE  35.  DEVELOPMENT  OF  LONGITUDINAL  TAYLOR  TIME  MICROSCALES 
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(Tig)  calculations.  Figure  36  is  presented  for  consideration.  Energy 
dissipation  rates  are  highest  in  the  recirculating  region.  They  drop 
to  essentially  zero  at  X/2Ro=1.00  and  then  gradually  increase  down¬ 
stream.  This  trend  provides  further  confirmation  of  the  fact  that, 
large  transfers  of  energy  from  the  mean  flow  to  the  turbulent  flow  are 
taking  place  at  high  rates  in  the  recirculation  region.  Finally, 
cylinder  rotation  at  u>=500  RPM  increases  the  rate  of  energy  transfer 
in  the  recirculating  region. 


IV  CONCLUSIONS 


The  primary  objective  of  this  investigation  is  to  determine  how 
well  the  photon  correlation/laser  velocimeter  system  performed  in  the 
seeded  flow  environment  of  a  two  dimensional  smoke  tunnel.  Several 
parametric  studies  were  conducted  to  evaluate  this  performance  and  the 
conclusions  resulting  from  them  are  listed  below. 

(1)  Marker  particle  concentration  ranges  used  in  this  investiga¬ 
tion  have  no  appreciable  effect  on  local  mean  velocities  or  turbulent 
intensities. 

(2)  Pinhole  diameters  must  be  correctly  chosen  for  the  flowfield 
under  study.  Small  pinhole  diameters  produce  incorrect  turbulent  in¬ 
tensity  data  in  regions  of  low  velocity  (e.g.  recirculation  regions 
behind  cylinders) .  Recall  that  this  is  a  consequence  of  the  apparent 
fewer  number  of  particles  traversing  the  control  volume  fringes.  In 
the  frees tream,  turbulent  intensity  data  with  varying  pinhole  diameters 
is  ambiguous  and  represents  an  incompatibility  between  the  smoke  in¬ 
jection  method  and  the  laser  velocimeter  system.  Further  study  of  this 
ambiguity  is  warranted.  Mean  velocity  measurements  are  not  affected  by 
pinhole  diameter. 

(3)  The  effects  of  frequency  shifting  in  recirculating  regions  or 
frees tream  regions  requires  further  study. 

(4)  Particle  lag  in  this  investigation  is  negligible. 

Conclusions  not  directly  related  to  the  primary  objective,  concern 

the  nature  of  the  circular  cylinder  flowfield  and  can  be  stated  as 
follows. 

(1)  Cylinder  rotation  effects  on  the  near  wake  (including  re¬ 
circulation  region)  are  confined  to  influencing  turbulent  intensities. 


Mean  velocity  effects  are  insignificant.  In  the  potential  flow  regions 
of  the  cylinder  flovfield  the  opposite  trend  is  noted.  Mean  velocities 
are  affected  and  turbulent  intensities  are  unaffected. 

(2)  The  integral  time  scales  are  larger  than  the  Taylor  time 
microscales  as  expected.  Both  time  scales  grow  slowly  downstream.  The 
fact  that  their  magnitudes  are  similar  indicates  that  the  energy  con¬ 
taining  region  and  dissipation  region  of  the  flow  energy  spectrum  are 
not  very  far  apart. 

(3)  High  rates  of  energy  transfer  from  the  mean  flow  to  the  turbu¬ 
lent  flow  are  taking  place  within  the  cylinder  recirculation  region. 
Observed  frequencies  of  the  integral  scale  and  Taylor  microscale  were 

on  the  order  of  100  KHz  to  200  KHz,  indicating  small  scale  turbulent 
structures  that  are  dissipative  of  energy. 

Based  upon  this  set  of  conclusions,  it  appears  that  the  laser 
velocimeter  system  and  the  two  dimensional  smoke  tunnel  are  highly  com¬ 
plementary  of  each  other,  but  further  work  is  required  in  the  inter¬ 
pretation  of  results  from  pinhole  diameter  effects  and  frequency  shift 
effects . 
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V  RECOMMENDATIONS 


The  following  recommendations  are  advanced  to  improve  the  quality 
of  investigations  performed  with  the  laser  velocime ter /smoke  tunnel 
system  and  to  serve  as  a  basis  for  future  investigations  with  the  same 
system. 

(1)  A  more  uniform  method  of  introducing  marker  particles  into 
the  flow  is  required  that  would  lower  freestream  turbulent  intensities 
and  make  measurements  in  the  freestream  independent  of  position  of  the 
laser  control  volume. 

(2)  A  traversing  table  should  be  designed  and  fabricated  that 
allows  the  laser  velocimeter  to  be  run  in  the  forward  scatter  mode  for 
this  flow  facility. 

(3)  The  electronic  data  processing  should  be  upgraded  to  eliminate 
the  time  consuming  hand  recording  of  data. 

(4)  Smaller  diameter  cylinders  should  be  tested  that  allow  measure¬ 
ments  to  be  taken  in  the  far  field  wake  for  comparison  with  theory  and 
prior  empirical  results. 

(5)  The  parametric  studies  begun  in  this  investigation  should  be 
expanded  to  include  further  values  of  Upg,  AF,  marker  particle  con¬ 
centration  and  co.  Pinhole  diameter  effects  with  the  existing  pinholes 
should  be  studied  in  greater  detail  at  different  flowfield  locations. 

(6)  The  scope  of  this  investigation  should  be  broadened  by  making 
measurements  of  power  spectral  densities,  intermittencies ,  convection 
velocities,  coherence  of  the  velocity  flowfield  and  pressure/velocity 
cross  correlations. 
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APPENDIX  A 


MEAN  VELOCITY  AND  TURBULENT  INTENSITY 
DATA  REDUCTION 

Longitudinal  mean  velocities  and  turbulent  intensities  are  computed 
from  a  typical  autocorrelation  function  shown  in  Figure  38.  This  auto¬ 
correlation  function  has  been  corrected  for  skewness  (Appendix  B) .  The 
computational  process  can  vary,  depending  upon  whether  a  frequency  shift 
has  been  applied  to  the  phase  modulator  or  not.  Recall  that  frequency 
shifts  are  required  in  regions  of  high  velocity  or  high  turbulence  and 
for  determining  flow  direction.  For  the  case  of  no  frequency  shifting 
(turbulent  intensities  generally  less  than  30%)  the  required  expressions 
are  as  follows  (Reference  12). 

For  the  mean  velocity  the  formula  is: 

U  =  S 

T 

X2 

where:  S  =  fringe  spacing  in  the  control  volume 

T  =  time  to  the  autocorrelation  function  second  maxima 
at  R2 

The  fringe  spacing  itself  is  calculated  according  to: 

s  =  l  X 

d  rj 

where:  A  =  distance  from  the  focusing  lens  to  the  control  volume 

d  *  separation  of  the  laser  beams  at  the  focusing  lens 
A  *  wavelength  of  the  laser  beam 
X]  *  refractive  index  (air  «  1.0) 

For  the  turbulent  intensity  calculations: 

„  ,  !  1/2  „ 

RMS  1  [-j  (r-1)  +  X  REF 

VTS  "  2,1  °TS 
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where:  r  =  I^-R^ 

VS 

R^=  autocorrelation  function  first  minima 

autocorrelation  function  second  maxima 

R^=  autocorrelation  function  second  minima 
r 

N  = _ £  =  number  of  fringes  in  the  control  volume 

S 

=  laser  beam  radius/f ringe  spacing 

In  this  investigation  the  term  in  the  expression  for  turbulent  intensity 

_3 

containing  N,  is  neglected  since  its  value  is  negligible  (1.65  x  10  ). 

In  the  case  where  a  frequency  shift  is  required  (high  velocity  or 

high  turbulence)  the  expressions  are  as  follows. 

For  mean  velocity: 

U  -  Fd  A 

2  sin  (2a) 

where:  Fd=Doppler  frequency  of  the  flow  without  a  frequency  shift 
A  =  wavelength  ot  the  laser  beam 

a  =  half  angle  between  the  laser  beams  at  the  control 
volume 

The  doppler  frequency  of  the  flow  without  a  frequency  shift  is  in  turn 
defined  as; 

Fd=Fd*  +  AF 
* 

where:  Fd  =  Doppler  freq^ncy  of  the  flow  with  a  frequency  shift 

AF  =  frequency  shift  applied  (sign  depends  on  flow  direc¬ 
tion) 

Then,  the  doppxer  frequency  of  the  flow  with  the  frequency  shift  is 
defined  as: 

Fd*  *  2U*  sin  (2a) 

A 
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where:  U 


mean  velocity  calculated  without  frequency  shift 


For  turbulent  intensity  calculations: 


U 

U 


RMS 

FS 


RMS 


REF 


U 

U 


REF 

FS 


where: 


RMS 

REF 


turbulent  intensity  calculated  without  frequency 
shif  t 


Finally,  for  the  formula  for  turbulent  intensities  to  be  applied  meaning¬ 


fully,  the  autocorrelation  function  must  be  corrected  for  skewness 


(Appendix  B) . 
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APPENDIX  B 


SKEWED  AUTOCORRELATION  FUNCTION 
CORRECTION  PROCEDURE 

If  the  measured  autocorrelation  function  is  skewed,  either  due  to 
a  small  number  of  control  volume  fringes  or  detection  of  background 
scattered  light,  it  must  first  be  corrected,  before  accurate  turbulent 
intensity  values  can  be  calculated.  Mean  velocity  calculations  are  not 
affected  by  the  autocorrelation  function  skewness. 

A  typical  skewed  autocorrelation  function  is  shown  in  Figure  37 
for  the  flowfield  location  X/2Ro=0.5,  Y/2Ro=0.0.  To  correct  for  skew- 
ness,  a  least  squares  fit  polynomial  regression  routine  (Reference  13) 
is  utilized.  A  second  degree  polynomial  is  fitted  to  the  first  three 
autocorrelation  function  maxima  and  another  second  degree  polynomial  is 
fitted  to  the  first  three  autocorrelation  function  minima.  The  result¬ 
ing  polynomial  coefficients  are  respectively  Ao,A^,A^  and  Bo*B^,B2  and 
the  polynomial  equations  of  the  form: 

F1(T)=Ao+A1T+A2T2 

F0(T)=B  +B, T+B.T2 
z  o  1  z 

The  arithmetic  mean  of  these  coefficients  is  taken  yielding 

C  =A  +B  ,  C  =A,+B,  ,  C  =A-+B„ 
o  O  Q  1  1  1  Z  l  Z 

2  2  2 
and  a  new  "mean  parabola"  of  the  form 

f3(t)*cd+c1t+c2t2 

is  constructed  as  in  Figure  37. 

This  curve  represents  the  new  reference  datum  for  the  skewness 
corrected  autocorrelation  function.  The  mean  curve  is  then  subtracted 
from  the  previous  maxima  and  minima  curve  fits  at  points  R^,  R^  and  R^. 
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The  resulting  differences  are  replotted  with  a  subsequent  shift  of  the 
mean  parabola  to  the  X  axis  (Figure  38) .  The  corrected  autocorrelation 
function  is  now  available  for  turbulent  intensity  calculations 
(Appendix  A) . 

The  failure  of  the  computational  schemes  presented  in  Appendix  A, 
when  applied  to  a  highly  skewed  autocorrelation  function,  are  a  con¬ 
sequence  of  the  fact  that,  in  these  cases  (R2“R3> ^^2*^1^ 9  resulting  in 
r  values  less  than  one  (undefined  turbulent  intensities).  This  tech¬ 
nique  is  applicable  to  all  autocorrelation  functions  above  a  skewness 
of  zero. 


P 


AUTOCORRELATION  FUNCTION  CORRECTED  FOR  SKEWNESS 


APPENDIX  C 


AMBIGUITY  NOISE  CORRECTION 


Ambiguity  noise  is  due  to  the  finite  size  of  the  laser  measuring 
volume.  The  resultant  finite  transit  times  of  marking  particles  as  they 
pass  through  the  fringes  in  a  random  manner  causes  an  apparent  turbulence 
to  be  seen. 

Reference  15  states  that  the  ambiguity  noise  from  the  autocorrelation 
function  at  T=0.0  should  be  a  constant  dependent  on  optical  setup  geometry. 

The  autocorrelation  function  is  determined  from  the  power  spectra  as 
derived  by  Lumley  (Reference  14) 

P'(F)=  [EXP{-(F+Fo)2a2/2}  (1) 

+EXP{-(F-Fo)202/2}] 

by  taking  its  Fourier  transform  yielding: 

R(T)*=EXP(-T2/2o2)COS(FoT)  (2) 

where:  Fo=  doppler  frequency 

a  *  a  constant 

The  ambiguity  noise  is  then  defined  as: 

AFo  «  __1_  (3) 

Fo  aFo 

where:  AFo  ■  root  mean  square  deviation  from  the  doppler  frequency 


According  to  Reference  15  an  expression  for  aFo  can  be  determined  from 
the  formula  for  the  doppler  frequency 

Fo  -  411  U  sin  a  (4) 

where:  U-  local  mean  velocity 

a*  half  angle  between  laser  beams  at  the  control  volume 
\m  laser  beam  wavelength 
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by  inserting  the  constant  a  yielding; 

aFo  =  AH  (aU)sin  a  (5) 

X 

From  the  expression  for  the  autocorrelation  function,  the  quantity  OU 

2 

can  be  related  to  the  distance  between  the  1/e  points  of  the  laser 
beam 

d*=2/2  (OU)  (6) 

which  is  also  equal  to  the  theoretical  diameter  of  the  laser  beams  at 
the  focal  point: 

d~  Xf  (7) 

r  JI 
o 

where:  f-  focal  length  of  the  lens 

dT=  distance  between  1/e  points 
rQ=  laser  beam  radius  entering  lens 

Equating  (6)  and  (7)  a  value  for  OU  can  be  determined,  which  is  then 
used  in  expression  (5)  to  find  aFo.  The  reciprocal  of  aFo  then 
represents  the  correction  due  to  ambiguity  noise  required  of  the  auto¬ 
correlation  function. 

For  this  investigation  the  value  of  AFo/Fo  was  determined  to  be 
3.9%.  This  value  is  then  subtracted  from  the  turbulent  intensities 
calculated  to  arrive  at  a  correction  for  ambiguity  noise.  The  actual 
steps  in  the  calculation  are  given  below: 

d  -2(6328xl0~10m)(l  m) 

1  11(0.0011  m) 

-  3.66  x  10"*4  m 
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APPENDIX  D 

INTEGRAL  TIME  SCALE  DETERMINATION 
The  integral  scale  is  defined  as  the  largest  turbulent  scale 
existing  in  a  flowfield.  In  the  photon  correlation  scheme,  the  auto- 
correlation  function  is  defined  in  terms  of  the  particle  transit  time 
between  the  fringes  of  the  control  volume.  For  the  integral  time  scale 
determination,  the  time  to  the  first  zero  crossing  of  the  autocorrelation 
function  is  used  (Figure  39»  Reference  8).  That  is, the  integral  velocity 
scale  is  first  determined  using  this  zero  crossing  time  as  follows: 


where:  S=  fringe  spacing 

T  ■  time  to  the  first  autocorrelation  zero  crossing 
IS 

This  represents  the  velocity  of  the  large  scale  disturbances  in  passing 
one  fringe  spacing.  From  this  velocity  and  the  fringe  spacing,  a  fre¬ 
quency  of  the  integral  scale  disturbances  can  be  determined 

Fis  m  UIS 


and,  proceeding  one  step  further  the  integral  time  can  be  determined  as: 


Integral  time  scales  computed  by  this  method  are  of  the  same  order 
of  magnitude,  but  larger  than,  the  Taylor  time  microscales  computed  by 
a  similar  method. 

Finally,  the  integral  time  scales  computed  by  this  method  should  be 
viewed  in  terms  of  trends  over  the  flowfield  and  not  in  terms  of  absolute 
values  due  to  the  approximate  nature  of  this  calculation. 
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FIGURE  39.  AUTOCORRELATION  FUNCTION  FOR  INTEGRAL  TIME  SCALE 

DETERMINATION 


APPENDIX  E 


TAYLOR  TIME  MICROSCALE  DETERMINATION 
The  Taylor  microscale  is  the  characteristic  turbulent  scale  present 
in  a  flowfield.  It  lies  between  the  integral  scale  and  the  Kolmogorov 
scale  Csmallest  scale  present  in  a  flowfield). 

The  first  step  in  the  determination  of  the  Taylor  time  microscale  is 
to  fit  a  parabola  to  the  autocorrelation  function  near  T=0.0  (Figure  40). 
Ideally,  this  parabola  should  correspond  closely  with  the  first  few  data 
points  of  the  autocorrelation  function  (Reference  16) . 

This  procedure  was  altered  for  this  investigation  by  fitting  a 
parabola  to  the  first  two  autocorrelation  function  points  excluding  the 
T-0.0  point.  This,  effectively  eliminates  any  inaccuracies  due  to 
ambiguity  noise  present  at  T=0.0. 

Due  to  the  low  number  of  autocorrelation  function  points  between 
T=0,0  and  T=T^  a  parabolic  fit  to  the  first  two  eligible  points  repre¬ 
sents  a  higher  percentage  of  the  autocorrelation  function  curve  itself 
then  just  the  accepted  first  few  data  points.  Since  the  Taylor  time 

microscales  are  based  upon  the  value  of  T  ,  it  can  be  expected  that 

SiS 

the  altered  procedure  will  result  in  Taylor  time  microscales  whose 
absolute  values  will  be  inaccurate.  However,  applying  the  procedure 
uniformly  to  all  autocorrelation  functions  should  result  in  observable 
trends . 

Once  the  parabolic  fit  has  been  accomplished  and  Tp  determined; 

Sis 

calculation  of  the  Taylor  time  microscales  proceeds  in  a  fashion  similar 
to  the  integral  time  scale  calculations.  The  Taylor  velocity  scale  can 


be  determined  as: 


s 


r 


where:  S  »  fringe  spacing 

T  =  time  to  the  intersection  of  the  parabola  with 
rlS  the  T  axis 

In  turn  a  Taylor  microscale  frequency  is  defined  as 


S 


and  the  desired  Taylor  time  microscale  is  the  inverse  of  the  Taylor 
microscale  frequency: 


As  with  the  integral  time  scale  calculations,  absolute  values  of 
T^  are  apt  to  be  subject  to  error.  Most  use  of  them  will  come  in  ob¬ 
serving  their  trends. 


L 
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APPENDIX  F 


PARTICLE  LAG  ESTIMATION 

The  question  of  whether  a  particle  in  a  flowfield  follows  the  fluid 
motion  or  not,  is  of  primary  concern  to  any  investigation  utilizing  a 
laser  velocimeter.  The  velocity  that  is  actually  measured  by  the  laser 
velocimeter  is  that  of  the  particles  and  not  the  fluid.  Regions  of  high 
acceleration  or  deceleration,  where  the  particle  may  not  follow  the  flow 
direction,  are  most  suspect  in  regards  to  this  question.  Hence,  flow 
around  a  two  dimensional  circular  cylinder  from  the  leading  edge  stagna¬ 
tion  point  to  the  upper  or  lower  midchord,  represents  the  most  critical 
area  of  concern  for  this  investigation,  since  the  highest  levels  of 
acceleration  are  achieved  here. 

An  analysis  of  particle  lag  (Reference  16)  begins  with  the 
equation  of  motion  of  a  particle  in  a  non-uniform  flow  under  the 
assumptions  that,  the  gas  density  is  much  less  than  the  particle 
density  and  body  forces  are  ignored  in  comparison  to  the  aerodynamic 
drag  of  the  particle: 

dUp  -  3  V  Re  CD  (U  -U  )  (1) 

dT  4  - V-  8  P 

PPD 

where:  p  *  coefficient  of  viscosity  (air) 

CD  -  CD(Re  ,  Mpg,  Knudsen  no.) 

P 

■  drag  coefficient  of  the  particle 

p  *  density  of  particle 
P 

D  -  particle  diameter 
11^  ■  gas  velocity 
Up  *  particle  velocity 
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# 


t  I 


Re  ■  Reynolds  number  based  on  particle  diameter 

P~  VVVP 

V 

From  Stokes  Law: 


e 


Substituting  equation  (2)  into  equation  (1)  yields: 


dt 


3  24 

4  R 


'S.  <yv 

ep  pp»2 


or 


dU 

_ E 

dt 


K  (U  -U  ) 
g  P 


where:  K  ®  18y  =  Stokes  coefficient 

V* 


(2) 


(3) 


(4) 


Particle  sizes  normally  used  for  flow  seeding  in  laser  velocimeter 
applications,  are  of  the  order  of  1  micron  diameter  or  less.  Particle 
sizes  resulting  from  vaporization  of  kerosene  are  assumed  to  be  larger 
and  for  this  particle  lag  estimation  a  particle  size  of  5  microns 
diameter  will  be  assumed* 

Computing  Stokes  coefficient  from  the  following  parameters: 
D*5ym®5xl0^m 
yflir  »  1*796  x  1(T5  kg/m  sec 

pp  -  1.4  x  103  kg/m3 


we  obtain: 

K  *  18(1.796  x  10~3kg/m  sec)  *  9.2  x  103 

(1.4  x  103kg/m3) (5  x  10_6m)2 


Assuming,  in  the  worst  case  situation  for  a  two  dimensional  cylinder. 
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that  the  flow  velocity  Increases  from  0  to  5.9  m/s ec  (U_,c  for  this 

investigation)  within  0.059  m(radlus  of  the  circular  cylinder),  then 

Ug  can  be  assumed  to  vary  linearly  as: 

U  -  3  X  where  0-100  (5) 

© 

If  the  particle  lag  can  be  expressed  as  U  -  U  -U  then  substitution 

—  g  p 

of  this  expression  into  equation  (4)  and  making  the  assumption  that 

U«U  results  in: 

-  g 

U  dU  U  dU  +  K  U  (6) 

g  _ g.  -  g  _=  - 

dx  dx 

Inserting  equation  (5)  into  equation  (6)  yields: 

02  X  -  0  X  dU  +  K  U  (7) 

dx 

Integrating  and  assuming  U  =  0  at  X  *  0,  then: 

U  -  J?_  X  *  0  Ug  (8) 

0+  K  0+  K 

Substituting  0  and  K  for  this  investigation  into  (8)  gives: 


U  -  100  U 


1.08  X  10-2  U 


100+9 . 2X10 J 

The  particle  lag  is  one  percent  of  the  gas  mean  velocity.  Since 
this  calculation  is  performed  for  the  region  of  highest  flow  accelera¬ 
tion,  the  particle  lag  in  less  accelerated  regions  can  be  assumed  to  be 
less  than  one  percent.  Therefore  particle  lag  is  considered  to  be 
negligible  for  this  investigation. 
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